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Abstract 



Static calculation and preliminary kinetic Monte Carlo simulation 
studies are undertaken for the nucleation and growth on a model sys- 
tem which follows a Prank-van der Merwe mechanism. In the present 
case, we consider the deposition of Ag on Au(lOO) and Au(lll) sur- 
faces. The interactions were calculated using the embedded atom 
model. The process of formation and growth of 2D Ag structures 
on Au(lOO) and Au(lll) is investigated and the influence of surface 
steps on this phenomenon is studied. Very different time scales are 
predicted for Ag diffusion on Au(lOO) and Au(lll), thus rendering 
very different regimes for the nucleation and growth of the related 2D 
phases. These observations are drawn from application of a model free 
of any adjustable parameter. 

Keywords: underpotential deposition, embedded atom method, 
dynamic Monte Carlo simulation, nucleation and growth, 2D phase 
formation. 
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1 Introduction 



The electrodeposition of a metal (Me) onto a foreign solid surface (S) is one 
of the most extensively studied subjects in surface electrochemistry. This 
electrochemical phase formation phenomenon is a key aspect in important 
technological processes such as electroplating and electrocatalysis. In the last 
years, the application of in-situ local probe microscopy (SPM) techniques 
such as scanning tunneling microscopy (STM) and atomic force microscopy 
(AFM) provided a powerful tool to observe the initial stages of these pro- 
cesses on an atomic level. In a number of cases, metal overlayers can be 
electrodeposited onto a foreing metal substrate at a potential that is less 
negative than the Nernst equilibrium potential of the 3D metal phase. This 
so-called underpotential deposition process (UPD) occurs in the undersatu- 
ration or underpotential range, given by convention as: 

AE = E- EsDMe > (1) 

where E is the actual electrode potential, E^^Me represents the equilibrium 
potential of the 3D metal phase, and AE is the underpotential shift. 

The UPD process has been well characterized for many systems with 
the SPM techniques and the preceding formation of metal phases of low 
dimensionality (OD, ID, 2D) have been analyzed [||, 0, H, ^, ^ |]. The stability 
ranges of these iD Me phases (« = 0, 1, 2, 3) can be formally described by 
Nernst-type equations: 

EiDMe = E^DMe (2) 

where aMe^+ denotes the activity of Me^"*" ions within the electrolyte, aiOMe is 
the activity of the iD Me phase, which is a constant for a condensed Me phase, 
and E^jjj^j^ represents the corresponding equilibrium potential. Usually, the 
electrochemical formation of the low dimensional condensed phases occurs 
under supersaturation conditions, whereas their dissolution takes place un- 
der undersaturation conditions. For an zD Me phase, the corresponding 
overpotential deposition (OPD) or supersaturation range, is given by: 

ViDMe = E — EiUMe < (3) 

where rjiDMe is the overpotential corresponding to the ?D Me phase. The 
formation of these iD Me phases is strongly influenced by surface defects 
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like kinks, vacancies, chemical impurities, monatomic steps, stacking faults, 
etc. The crystal surface can also be considered as a 2D crystal imperfection 
and plays an important role in the UPD phenomena. Thus, the stepwise 
formation of low-dimensional Me phases on a substrate in the underpotential 
range are generally characterized by the formation of OD and ID phases at 
relatively high underpotentials, followed by the formation of expanded and 
condensed 2D phases on atomically fiat terraces at lower underpotentials. 

In a previous work [0, we considered a theoretical model suitable to 
analyze the existence of these different phases when a metal is deposited 
on a foreign substrate with defects. We employed a lattice model along 
with a grand Canonical Monte Carlo simulation to study the thermodynamic 
features of the deposition process. One of the novel features of this model 
consisted in the realistic modelling of the interactions between the particles 
of the system, which was made using the embedded atom method (EAM). In 
this way, we introduced in a simple lattice model the many body features of 
the metallic binding which are indispensable for a reliable description of the 
metal-metal interactions. In agreement with the expectation drawn from the 
experience with experimental systems, we observed the correct sequence of 
appeareance of 0D,1D, and 2D phases when changing the chemical potential 
of the system in a direction corresponding to polarization of the surface 
towards negative potentials. Furthermore, the different regions of existence 
of the phases were separated by chemical potentials of a few hundreds of meV, 
whose values are of the order of magnitud of those observed experimentally. 

In some cases, the 2D Me phase formation is not reached after sequence 
of appearance of OD and ID phases, but proceeds by the formation of 2D 
expanded commensurate Me adlayers which are transformed at lower under- 
saturation into condensed metal monolayers via a first order phase transition. 
Such phase transitions are observed, for example, during the Ag UPD on both 
Au(lll) and Au(lOO) surfaces |, |]. 

In summary, UPD systems are very complex ones, as they involve phases 
of different dimensionalities which are mainly characterized by: 

1. Specific adsorption of metal ions 

2. Electronic charge transfer 

3. Metal adatom-adatom and substrate-adatom interaction, which can- 
not be rigorously disentangled due to the non-additive property of the 
metal-metal interaction. 
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4. Undersaturation and/or supersaturation conditions with respect to the 
several {iD) phases that may appear on the system. 

5. Crystallographic misfit, responsible in many cases for the lack of reg- 
istry between the adsorbate and the substrate. 

6. Cosorption and competitive sorption phenomena of anions as well ox- 
idic species leading also to the formation of low-dimensional anion or 
oxidic systems in the corresponding undersaturation ranges. 

7. The influence of the zero charge potential of the substrate on the electric 
field conditions at the interface and on adsorption and phase formation 
phenomena. 

8. Surface site exchange processes leading to the formation of surface 
Me-S alloys. 

In principle, a general description of the UPD phenomena would require 
a model incorporating all these items, something that is possible in principle 
with the current status of the developing of computers but still a formidable 
task for the small community of theoretical electrochemistry. However, more 
simple models may be formulated, including some of the points mentioned 
above, and a number of properties may be inferred for those systems whose 
salient features correspond to the model selected. 

The key features of the present work concentrate on a reliable descrip- 
tion of the metal-metal interaction by the EAM, and a lattice model which 
allows the consideration of relative large systems. Thus we can make the 
following comments in relation with the points mentioned above: Points 1), 
2), 6) and 7) cannot be directly addressed in the present model, since this 
would require a detailed description of the metal/solution interface, that is 
only possible using first-principles calculations. Therefore, we shall only deal 
with the deposition of atoms at different rates, without being concerned on 
how these rates are related to the corresponding potential changes at the 
interface. Point 3) is one of the strenghts of the present formulation, since 
the interactions in the EAM have been devised to account for the properties 
of the metals involved and their alloys. With respect to point 4), in the 
present work we concentrate on the deposition process, neglecting adatom 
desorption. So we set supersaturation conditions for the formation of the 
several low dimensional Me phases. Point 5) is a very complicated one, since 
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the crystallographic misfit induces in many cases the existence of incommen- 
surate structures, which cannot be addressed by lattice models. Thus, the 
way to overcome this problem is by choosing a system with a negligible crys- 
tallographic misfit. In this respect, the deposition of Ag on Au appears as 
an optimal candidate. Concerning point 8), site exchange with the substrate 
could be considered within the present model, but preliminary studies by 
molecular dynamics show that this effect is not of primary importance in the 
submonolayer range. 

According to the previous discussion, we perform in this work model cal- 
culations for Ag deposition on Au(lll) and Au(lOO) surfaces, modelling the 
metal-metal interactions by the embedded atom method. These systems ap- 
pear as optimal candidates due to a negligible crystallographic misfit. On 
the other hand, both Ag and Au atoms are mobile in our model, so that 2D 
alloys could naturally appear during the course of the simulation but this 
2D alloying process has not been observed. Therefore, the aim of this work 
is to set the basis for a systematic study of a system involving electrochem- 
ical phase formation in a Frank-van der Merwe system and emphasize the 
role of metal-metal interactions in the surface processes accompanying this 
phenomenon. 

2 Some experimental facts of the system 

Au(/i/c/)/Ag+ 

The system Au(/ifc/)/Ag"^ is a typical example for Me UPD on a foreign sub- 
strate S, with strong Me-S interaction but negligible Me-S misfit {cIq au = 
0, 2884 nm, d^^Ag = 0, 2889 nm). Thus, the Frank- van der Merwe or layer by 
layer growth mechanism is expected to operate in this system. Our previous 
studies p, ^ have demonstrated that the Ag UPD on Au(lll) and Au(lOO) 
occurs in the underpotential range < AE/voN < 720. For the system 
Au(100)/Ag+, the formation of an expanded Au(lOO) -c(v^ x 5y2)R45°Ag 
overlayer is observed within the potential range 200 < AE/mV < 550. At 
lower underpotentials a limited growth of steps occurs and 2D Ag islands 
are formed involving simultaneous 2D nucleation and growth at steps and 
fiat terraces (Figure 1). The atomic structure on top of islands as well as 
terraces is quadratic with an interatomic distance of do,Ag = 0.29 ± 0.01 nm. 
These morphological and atomic observations indicate that the expanded 
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structure transforms into a Au(lOO) — (1 x l)Ag phase via a first order phase 
transition. The steps of the 2D Ag islands grow slowly depending on the 
potential, until a complete Ag UPD monolayer is formed at AE < 15 mV. 
On the other hand, in the system Au(lll)/Ag"'" the experiments showed an 
expanded Au(lll) — (4 x 4)Ag structure in the range 50 < AE/mV < 500 
which is transformed to a Au(lll) — (1 x l)Ag phase at lower underpoten- 
tials. The interatomic distance observed {do^Ag = 0.28 ± 0.01 nm) indicates 
the formation of an hexagonal close-packed Ag monolayer. It is important 
to note that, in this case, the formation of the condensed 2D Ag phase on 
Au(lll) occurs preferentially at monatomic steps of the substrate, without 
the formation of 2D islands (Figure 2). The atomic structures and mor- 
phologies described above were observed in both sulphate and perchlorate 
solutions. In addition, electrosorption valency measurements have indicated 
that coadsorption or competitive adsorption of sulphate or perchlorate an- 
ions can be excluded. Nevertheless, as was suggested p, , a nearly constant 
anion coverage in the entire Ag UPD range cannot be excluded because Ag 
UPD on An{hkl) occurs at positive potentials with respect to the potential 
of zero charge of Au{hkl) and Ag{hkl). 

From this brief experimental outlook, it is clear that the model we are 
presenting in this work is still missing an important number of features to 
reflect the experimental situation. However, as we shall see below several 
aspects of the experimental system coincide with the predictions of the model. 



3 Model and simulation method 



3.1 Interatomic potential 

The choice of proper interatomic potentials to perform the simulations is one 
of the key elements of the model. Several methods are available to calculate 
the total energy of a many-particle metallic system, with a computational 
effort comparable to that of a pair potential [10|. It is worth mentioning the 
embedded atom method (EAM) |Tl|, the N-body potentials of Finnis and 
Sinclair 0], the second-moment approximation or Tight-binding (TB) ITT 



and the glue model (GM) |]T3|. In this work we employ the embedded atom 



method |TT[ because it reproduces important characteristics of the metallic 
binding that cannot be obtained using simple pair potentials. 

The EAM considers that the total energy Utot of an arrangement of 
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particles may be calculated as the sum of energies Ui corresponding to indi- 
vidual particles 

N 

Utot = Y.U^ (4) 
1=1 

where Ui is given by 

Fi is denominated embedding function and represents the energy necessary 
to embed atom i in the electronic density ph,i at the site at which this atom 
is located, ph^i is calculated as the superposition of the individual electronic 
densities Piirij) : 

Ph,i = Epi(rij) (6) 

Thus, the attractive contribution to the EAM potential is given by the 
embedding energy, which accounts for many-body effects. On the other 
hand, the repulsion between ion cores is represented through a pair potential 
Vij(rjj), which depends only on the distance between the cores rj^: 

^ Zi{rij)Zj{rij) 

where Zi{rij) may be envisaged as a sort of effective charge, dependent on 
the nature of the particle i. 

The EAM has been parametrized to fit experimental data like elastic 
constants, dissolution enthalpies of binary alloys, bulk lattice constants and 
sublimation heaths |Tl|. Pair functionals have been successfully used for 



surface diffusion studies and adsorption of metals on metallic surfaces ||I4 

m ra. 



3.2 Lattice model 

Lattice models for computer simulations have been widely used in studies of 
nucleation and growth, because they allow simulations with a large number 
of particles at a relatively low computational cost. The reason for this advan- 
tage is the use of fixed rigid lattices that restrict enormously the number of 
possible configurations for the adsorbate as compared with a model where in 
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principle all positions in space are available. While continuum Hamiltonians 
are much more realistic in those cases where epitaxial growth of an adsor- 
bate leads to incommensurate adsorbed phases Jl^ or to adsorbates with 
large coincidence cells, the use of lattice models may be justified on the basis 
of experimental evidence or continuum computer simulations that predict a 
proper fixed lattice geometry. In the present case, we have strong evidence 
from simulations within the canonical Monte Carlo method |jl8| which in- 
dicates that at least one of the phases occurring during Ag underpotential 
deposition on Au(lll) and Au(lOO) possess a pseudomorphic structure. In 
fact, our continuum MC simulations showed that Ag monolayers adsorbed 
on Au(lll) and Au(lOO) spontaneously acquired a (1 x 1) coincidence cell 
in agreement with the experimental finding at low underpotentials. For this 
reason, we shall employ here a lattice model to represent the square (100) 
and the hexagonal (111) surface lattices in kinetic Monte Carlo simulations. 
Square lattices and hexagonal lattices of different sizes with periodical bound- 
ary conditions are used in the present work to represent the surface of the 
electrode. Each lattice node represents an adsorption site for a Ag or a Au 
atom. 

Two further approximations in the present model must also be mentioned. 
First, we neglect the effect of the presence of solvent molecules. This approx- 
imation should not be critical as long as the partial charge on the adatoms 
is small, thus minimizing the ion-dipole interactions. Second, we also ne- 
glect effects related to specific anion adsorption. This may not be true for 
the experimental system, but thin layer twin electrode experiments gave no 
indication of any change in the amount of anion adsorption upon building 
the adsorbed monolayer 0. On the other hand, the possibility of a nearly 
constant anion layer on the bare substrate Au(lOO) as well as on a Ag UPD 
modified Au(lOO) has been indicated 0. In the case of other systems, the 
adatom-anion interactions have been shown to be very important |]19| in 



playing a decisive role for determining the energetics of the system. This has 



been recently analyzed by some of us in thermodynamic terms |20 



3.3 Calculation of Gibbs energy of cluster formation. 

The initial stages of a phase formation are of primary importance, because 
the competition of different processes occurring during this phenomenon de- 
termine in many practical cases the final structure of the deposit. The ther- 
modynamic foundations of the initial stages of bulk phase formation are well 
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settled and have been thoroughly discussed in advanced books on the field 
m. For this reason, we shall address here briefiy the topics relevant for 
the present calculations. The free energy change AG{N) produced when 
a three-dimensional cluster of N atoms is formed on a surface at a certain 
overpotential rjsDMe can be written as: 

AG{N) = -Nze \r]sDMe\ + $(iV) (8) 

where z is the valence of the deposited ions, e is the elementary charge 
and $(A^) is the energy consumed for the formation of the new interface 
boundaries. In the 3D case, this quantity is related to the specific surface 
energies cTj of the facets of the cluster, the specific interface energy aj* and 
the surface energy of the substrate according to: 

HN) = (^iA + Aj* {(7j* - asub) (9) 

where At denotes the respective surface areas. ^{N) can also be calculated 
within the atomistic model, as long suitable potentials for the system are 
available. If we called E^+s the energy of a system consisting of a cluster 
of atoms of the metal Me on a substrate S, Es the energy of the free 
substrate, and Ei, the binding energy per atom of the bulk adsorbate metal, 
$(A^) would be given by: 

^N) = E,^s -Es-N E, (10) 

If the energy of the atoms in the cluster was equal to that in the bulk, Eq.p!0| 
would yield ^{N) = 0, and a vanishing small overpotential would be enough 
for the growth of the new phase according to Eq. p. 

In case of 2D nucleation, which is the current subject of our interest, an 
equation similar to is valid. However, instead of defining the overpotential 
with respect to the deposition of bulk Me, it is more suitable to refer to it as 
the thermodynamic deposition potential of the monolayer. Thus, the excess 
energy $2d(^) "will be calculated from: 

^MN) = E^^s -Es-N E^^ (11) 

where E^^^ is the binding energy of the Me atoms in the monolayer. In 
consequence, $2d(-^) takes into account the excess energy connected with 
the occurrence of the border of the 2D islands. 
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For a given number of atoms A^, several values of ^2d{N) can be found 
depending on the geometry of the 2D cluster. For each we have only 
considered the energy corresponding to most stable configuration. This was 
obtained by means of a simulated annealing procedure. This technique has 
often been used to obtain minimal energy structures or to solve ergodicity 
problems. A suitable way to implement this technique is through the canon- 
ical Monte Carlo method at different temperatures. In other words, a given 
number of atoms, A^, is selected, and a simulation allowing the motion of the 
atoms on the surface is started at a high initial temperature To, in the order 
of 10^ K. The system is later cooled down following a logarithmic law: 

Tf = ToK^^y^'-^ (12) 

where Tj is the final temperature, Ncydes is the number of cooling steps 
and i^' is a positive constant lower than 1. A few hundreds of thousands of 
Monte Carlo steps are run at each temperature in order to allow an extensive 
exploration of the configuration space and the simulation stops when Tf is 
reached. 



3.4 Model for Dynamic Monte Carlo simulations 

Although Monte Carlo methods have been traditionally related to the study 
of equilibrium properties, it is possible, if some conditions are fulfilled, to 
use them to compute the time evolution of a system. The foundations of 
dynamical Monte Carlo simulations have been discussed by Fichtorn and 
Weinberg |^ in terms of the theory of Poisson processes. According to this 



method, the sampling of the system must implicate transition probabilities 
based on a reasonable dynamic model of the physical phenomena involved. 
Besides fulfilling the detailed balance criterion, the transition probabilities 
should reflect a "dynamic hierarchy " related to the processes taking place 
in the system. 

In this preliminary study we are mainly interested in nucleation and 
growth phenomena, therefore, we shall consider only two types of processes 
related to the growth of the new phase: adsorption of an adatom on the 
surface and its diffusion in different environments. Since we are neglecting 
the desorption processes, the results of our model will be valid for relatively 
large overpotentials for the deposition of the Ag atoms. 

To illustrate the model used for the dynamic calculations, we show in 
Figure 3 a-c the diffusion of a Ag atom on a Au(lOO) surface in three different 
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environments. In Figure 3a, a Ag adatom jumps between two equivalent sites 
without any near neighbors. In Figure 3b, it has one Ag neighboring atom 
in the initial position and none in the final state and in figure 3c the jump is 
again between two equivalent positions with one Au nearest neighbor. The 
corresponding potential energy curves calculated by the EAM at each stage 
of the diffusion are shown in Figure 3d. From these curves we obtain the 
two important quantities required to perform the dynamic MC simulation: 
the attempt frequency for overcoming the diffusion barrier, calculated for 
the curvature of the potential energy surface at the initial state, and the 
activation energy, by taking the difference between the energy at the saddle 
point and the corresponding value at the initial minimum. These curves were 
constructed for all possible environments and involved 729 and 6561 energy 
curves in the case of the (100) and (111) faces respectively. 

In the case of the adsorption rates, we assumed the same rate for all free 
sites. A more realistic calculation should also take into account different 
rates on different environments. This would require a complete knowledge 
of the degree of solvation of the discharging ion and detailed information of 
the different transition states in different environments. This information 
is not available for the present system but could be easily introduced in 
our modelling. Experimental evidence show that electrodeposition of metal 
ions on native or foreign substrates take place preferentially at OD and ID 
surface inhomogeneities such as kink and monoatomic steps. In the case 
of foreign substrates, OD metal clusters are preferentially formed at these 
inhomogeneities. We have addressed the role of different types of defects by 
means of thermodynamic Monte Carlo studies in previous work, confirming 
this expectation. In the present work, we study the deposition phenomena 
from a kinetic point of view, analyzing the influence of the deposition rate 
on the formation and growth of the 2-D phase on a foreign surface with a 
relatively small number of defects. 

4 Results and discussion 
4.1 Static calculations 

4.1.1 2D Ag cluster formation on Au(lOO) and Au(lll) surfaces. 

Figure 4 shows the Gibbs energy of 2D Ag cluster formation on Au(lOO) as 
a function of size N, calculated according to the method described in section 
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3.3. Since we are using an atomistic model, our curves are discontinuous. For 
this reason, we have fitted our curve for rj2DAg = by means of a square root 
law, and later employed it to draw the corresponding continuous lines for 
\V2DAg\ > 0. With this information we were able to calculate the size of the 
2D critical clusters for different overvoltages, as reported in Table la. The 
sizes of the critical clusters are relatively large as compared with the values 
estimated for nucleation of Ag on a Ag(lOO) quasi perfect surface |jl|, given in 
Table 2b. In other words, the overpotentials for nucleation are considerably 
larger when Ag nucleates on the foreign Au(lOO) surface. 

Figure 5 shows the Gibbs energy of a 2D Ag cluster formation on Au(lll) 
at different overpotentials. It can be noticed that for a given overpotential, 
the maximum in the curve appears at larger N than in the case of the Au(lOO) 
surface, that is, the energy for the formation of border atoms from atoms 
located in the monolayer is larger for the adsorbate on the (111) face. This 
fact can be understood in terms of the relative energy change for this process. 
In the case of an adsorbate in the (100) monolayer, it has eight nearest 
neighbors: 4 Au substrate atoms plus 4 Ag atoms belonging to the monolayer. 
This coordination changes to seven (4 Au + 3 Ag) for a border atom. In the 
case of the (111) adsorbate, the change is from nine (3 Au + 6 Ag) to seven 
(3 Au + 4 Ag) at the border, that is, the change of coordination to produce 
a border atom from the monolayer is more important in the case of the more 
compact adsorbate, originating a stronger energy change. This is somehow 
the 2D analogous of the qualitative justification why more open faces exhibit 
a larger surface energy than the more compact ones. 

4.1.2 Diffusion of single Ag atoms on Au(lOO) and Au(lll) 

The growth of a 2D phase after nucleation involves a number of elementary 
steps that should be considered in a general formulation of the phenomenon. 
This has been analyzed in detail in specialized texts on the field[|l|, hence we 
simply enumerate them here: 

1. Bulk diffusion of the discharging ions. 

2. Adatom electrodeposition and electrodesorption. 

3. Surface diffusion of the discharged adatoms and diffusion of adatoms 
along the steps. 
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4. Direct attachment to steps or kink sites. 



The relative rates of these processes govern the overall rate and it is of 
key importance for the modelling of each system to determine the one(s) 
ruling the whole process. For example, in the case of the growth of a 2D 
Ag cluster on Ag(lOO) quasi perfect faces, step 2 is very fast under the usual 
deposition conditions yielding an unusually high exchange current density of 
Ag atoms. In this way, the average random surface displacement of adatoms 
during their stay at the surface is small and surface diffusion plays only a 
subordinate role. 

The high exchange current density of Ag atoms has also dramatic effects 
during the deposition of the second Ag monolayer on Au(lll), producing very 
noisy STM pictures under these conditions. On the contrary, images obtained 
for the growth of Ag islands on Au(lOO) are stable and clear 0] (Figure lb, 
Ic), denoting that the exchange current for Ag atoms on Au surfaces is 
considerably lower. On these grounds, we assume in our simulations and the 
following considerations that atom deposition occurs under conditions where 
atom exchange is negligible. 

The diffusion coefficients Dhki of a Ag single atom on the Au surfaces can 
be calculated from: 



where a is the distance between adsorption sites, u is the attempt frequency, 
Ea is the activation energy for diffusion and n = 4 or 3 for the (100) and 
(111) faces respectively. We obtained Diqq = 3.3 x lO^A^s'^ and Dm = 
1.2 X 10^°A^s~^, indicating that the diffusion on the (111) face is much faster 
than the same process on the (100) face. These figures already give us a 
very important hint to predict the type of nucleation and growth behavior 
in these two faces under the usual experimental conditions in voltammetry. 
In fact, let us consider an experiment where the reduction current density 
producing adatoms, i, is of the order of 1 fiA/cm'^, and let us assume that 
the surface presents terraces of the order of 40 x 40 nm^. The number of 
atoms originated at the terrace per unit time will thus be of the order of 10^ 
. We will denote this number with -. Thus, once an atom is deposited, 
it has the chance to diffuse the length: 




(14) 
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before any other atom is deposited on the average on the same terrace. If 
we replace in this equation the Dh^i obtained, we get diQQ ^ 8 A and dm ^ 
1.5 X lO'^A respectively. Thus, while on the Au(lOO) face the deposited Ag 
atoms may meet other atoms for the nucleation and growth of the 2D phase, 
the Ag atoms deposited on the Au(lll) surface will diffuse unhindered to 
the border of the terraces and produce the 2D growth there. Since the 
diffusion barriers are usually higher on (100) surfaces than on (111) ones, we 
expect that this should be a rather general result, also valid for other upd 
systems. It is worth mentioning that, as shown in section 2, the experimental 
Ag/Au(hkl) upd system presents a similar behavior, in that the growth of 
the 2D phase occurs at the border of steps on the (111) face, and in the form 
of islands on the (100) face, besides a hmited growth at steps. 

4.2 Dynamic Monte Carlo simulations. 
4.2.1 Ag on Au(lOO) 

We discuss here the results obtained with the simulation 2D box shown in 
Figure 6, which is relatively small as compared with real terraces, but should 
yield the correct qualitative features. We consider 1600 adsorption sites 
distributed over a 115 x 115 square lattice, delimited by two rows of 
monoatomic high Au steps on each side. Periodic boundary conditions were 
applied in the direction parallel to the borders. Real surfaces present a 
more complex topology, with different kinds of OD and ID inhomogeneities 
like kinks or roughened steps. However, our simulation could also represent a 
situation where kink or steps are already decorated by OD or ID silver phases. 
We only allowed for the deposition of Ag atoms on the terraces at different 
rates Vads- We varied Vads between 10"^ s'^ and 10^ per free adsorption 
site, corresponding to potentiostatic conditions with an initial current density 
(© 0) ranging between 2x10"^ ^A/cm? and 2 x lOV^/cm^- Although 
we made several sets of simulations, we discuss here some representative 
examples to obtain a physical insight into the phenomena taking place in 
each case. 

For the slowest rate (10~^ s~^) we typically observed initially a single par- 
ticle diffusing on the surface towards the Au border step and then remaining 
attached to the border and diffusing parallel to it. Diffusion along the borders 
is relatively fast with respect to surface diffusion, since Dl^^^^ = 1.7 x 10^ 
A^s~^. About 15 seconds later, a second particle entered and diffused to the 
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opposite side, followed by others that behaved in the same way. Two Ag 
stripes parallel to the Au borders grew in this way, being completed at about 
650 seconds. After that a second row started to grow. Diffusion of Ag on the 
Ag border also showed a diffusion coefficient of the order of 2 x 10^ A^s~^. 
In this simulation, no Ag island formation was evident on the surface and 
the growth occurred at the Au border. The third Ag row was completed at 
about 1800 seconds. 

When the adsorption rate was Vads = 10^'^ s^^, all the particles deposited 
up to 50 seconds were found to diffuse to the borders. At that time, a dimer 
was formed, which initially diffused on the surface and later was fixed by 
addition of more particles. At 170 s a trimer was formed, giving place to a 
second island. At 700 seconds, the coalescence started to take place between 
the growing islands and the growth at the Au borders. At 7250 seconds, only 
a single vacancy remained on the surface, diffusing with Dyacancy = 4.2 x 10^ 

The run with Vads — 10~^ presented several interesting features, since 
at this deposition rate (corresponding to ^ 2 fiA/cm'^) the number of islands 
appearing in our simulation box allows for statistical analysis. For example, 
the number of islands was followed as a function of time. In the simulation 
probed, the number of islands reached a maximum of 7 at about 15 seconds 
(Figure 7), decreasing later and generating islands of vacancies ai t ^ 400 s. 
Att = 630 s,only a vacancy remained diffusing on the surface and the surface 
was completely covered by Ag atoms at t = 860 s. In Figure 8 we show 
snapshots of a simulation with this adsorption rate. 

In the case of deposition rates of 10""^ and larger, a coexistence of 
growing islands and diffusing particles is observed from the very beginning 
of the simulation. Typical times for the completion of the monolayer are re- 
ported in Table 2. It can be noticed that they follow the expected logarithmic 
law. 

Figure 9 shows the number of islands as a function of the coverage degree 
for different deposition rates. Each of the curves represents an average over 
six simulation runs. It can be noticed there how the number of islands 
increases with the deposition rate, with a maximum in the range 0.1 < G < 
0.2 which is progressively better defined at the higher rates. 

Unfortunately there are no experiments available for the present simu- 
lation conditions, consequently a direct comparison is not possible. STM- 
voltammetric experiments where the deposition rate is of the same order of 
magnitude as some of our results (i ft; 2 /X74/cm^, corresponding to Vads — 
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lO^^s-i)), show the occurrence of growing islands generated at times of 
the order of a few seconds, which correspond to the simulation times where 
the islands appear in our simulations. It must be stressed, however, that the 
deposition history in the experiment is completely different, since the pseu- 
domorphic phase appears after the voltammetric formation of an expanded 
Au(lOO) - c(y2 X 5x/2)i?45°Ag ph ase which is not considered in the present 
model. 

4.2.2 Ag on Au(lll) 

The simulation box was analogous to that of the Au(lll) face, and the 
adsorption rates were varied between 10^^ < Vads/s~^ < 10^. Since the 
activation barriers on this surface are very low, the absolute-rates model 
employed to calculate the transitions of atoms on the surface should deliver 
only qualitative predictions. Furthermore, some transitions on these surfaces 
were not activated at all. We assigned to these processes an arbitrarily high 
rate, so they made a negligible contribution to the time accumulation in the 
simulation and we were able to follow the relatively slow processes occurring 
on the surface. 

In the case of the Au(lll) surface, the adsorption sites constitute an 
hexagonal lattice. Two kinds of adsorption sites occur, usually denoted as 
'hep' and 'fee', depending on the position of the adsorbing atom with re- 
spect to the second lattice plane of substrate. In our model these sites have 
practically the same energy, but the lattices they define cannot be filled 
simultaneously due to steric hindrances. Thus, under some simulation con- 
ditions two domains (one fee and one hep) of the adsorbate appear, with a 
corresponding domain border. 

For Vads ^ 10^ s~^, all atoms reached the border of the box. However, 
since we are dealing now with the (111) face, the borders are no longer 
equivalent. One of them exhibits a facet with square symmetry (that we 
shall label {100}) and the other presents three fold sites {111}. In the case of 
the {100} border, we obtained straight Ag stripes parallel to the step, while 
a triangle like decoration was found in the case of the {111} step (Figure 
10). This different type of growth occurs in our model on kinetic ground 
as a consequence of the fact that only jumps between nearest neighbors are 
allowed. Molecular dynamic calculations would be very helpful to elucidate 
this point. 

For Vads = 10^ we observed more than one Ag atom diffusing simul- 
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taneously on the surface, building in some cases dimers that also diffuse. 

For Vads = 10^ we observed the formation of islands on the surface, 
initially triangular like. Domains of sites corresponding to adsorption on fee 
or hep sites are formed. The surface becomes completely filled at t = 7x 10~^s 

For 10^ < Vads/s—1 < 10^ hep and fee domains are formed on the surface 
(Figure 11), with the domain border fluctuating rapidly within one near 
neighbor distance. 

According to the simulation results for this phase, we conclude that for the 
deposition rates under which the experiments are undertaken (corresponding 
to 0.1 < Vads/s—1 < 1), no islands should be formed on the surface. Sur- 
face diffusion is fast and allows the atoms to reach neighboring steps on the 
surface, where nucleation and growth occur. Although the Ag deposition on 
Au(lll) seems to be a rather complex process, including the formation of a 
number of expanded faces, we think that these findings, that are based on the 
fast diffusional behavior of Ag on Au(lll), should be at least qualitatively 
valid. 

5 Conclusions 

Inspired in puzzling experimental results, we have modeled the underpoten- 
tial deposition for a Frank-van der Merwe system on (111) and (100) faces by 
means of a lattice model employing a realistic potential for the metal-metal 
interactions. From a thermodynamic viewpoint, the present static calcula- 
tions indicate that the overpotentials for Ag nucleation on Au will be larger 
for the (111) than for the (100) face. The reason for this behaviour is the 
larger energy that is required to generate the borders of the growing Ag 
clusters in the former case. 

If a negligible exchange current for adatoms is assumed, kinetic simula- 
tions by means of dynamic Monte Carlo indicate that the nucleation and 
growth should take place with characteristic time of tenths of seconds in the 
case of the Au(lOO) face and of the order of milliseconds in the case of the 
Au(lll) face. Therefore, if Ag is deposited at intermediate rates on both sur- 
faces, island growth should occur on Au(lOO) but not on Au(lll). It must 
be stressed that this behaviour is a purely kinetic effect. On thermodynamic 
grounds, the growth of the 2D phase should always start at steps and not on 
the terraces. 

It must be emphasized that these conclusions are all drawn from a model 
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which is free from any adjustable parameters. A desirable improvement of 
the present model is to consider the electrodesorption of adatoms in different 
environments that involves a model for electron transfer in this system. Work 
in this direction is in progress. 
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7 Tables 



Table 1: 

a) Size Ncru of the 2-D critical cluster of Ag on Au(lOO) for different 
overvoltages r/, calculated from the embedded atom method and a simulated 
annealing procedure. 





\ri2DMe\ /eV 


313 


10 


78 


20 


34 


30 


20 


40 


9 


60 


5 


80 


3 


100 


1 


170 



b) Sizes of critical clusters estimated for nucleation of Ag on a Ag(lOO) 
quasi perfect surface , as reported in ref. ||l| 



Ncrit 




V2DMe 


/eV 


67 
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25 


10 
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Table 2: Average times required for the completion of the monolayer 
Tmon ior different adsorption rates per site Vads- These times correspond to 
values averaged over 6 simulation runs. 





^ niun / 




6.7x10'! 




8.3x10^ 


10-2 


7.50 xl02 


10-1 


70 


1 
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101 


0.69 


10^ 


0.067 
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8 Figure Captions 



Figure 1 In situ STM images of Ag UPD in tlie system Au(100)/5 x lO^^M 
AgC104 + 5 X HCIO4 at T = 298K. a) = 700 mV, b) = 200 

mV, c) AE = 20 mV. lup = 20 nA. 

Figure 2 In situ STM images of Ag UPD in the system Au(lll)/5 x 
lO-^M AgC104 + 5 X 10-^M Ha04 at T = 298K. a) AE = 700 mV, b) 
AE = 278 mV, c) AE = 27 mV. lup = 15 nA. 

Figure 3: a-c) Sample environments for the motion of a Ag atom on a 
Au(lOO) surface d) potential energy as a function of the distance along the 
diffusion path for the environments shown in a-c. 

Figure 4 Gibbs energy of 2D Ag cluster formation on Au(lOO) as a 
function of size N at different ovcrpotentials r]2DMe- 

Figure 5 Gibbs energy of 2D Ag cluster formation on Au(lll) as a 
function of size N at different overpotentials r]2DMe- 

Figure 6: One half of the simulation box employed to represent Ag 
nucleation and growth on a Au(lOO) surface. 

Figure 7: Number of islands as a function of time for a deposition rate 
of lO-^s-^ 

Figure 8. Snapshots of a simulation with the Au(lOO) surface, Vads = 
IQ-^s'^. ti=1.7 s, t2=5.7 s , t3= 19 s, t4= 42 s, tg^ 76 s, te^lOl s. Grey 

squares are Ag adatoms, black squares represent Au border atoms. 

Figure 9 Number of islands as a function of the coverage degree for 
different deposition rates v. Each curve represents an average over six simu- 
lation runs. 

Figure 10. Snapshots of a simulation with the Au(lll) surface Vads — 
10^s~^. Simulation times arc ti= 0.8 s, t2=4 s, t3= 6.6 s , t4= 13 s , t5= 
42 s, t6=87 s. The upper border corresponds to the {100}step and the lower 
border to the {111} step. 

Figure 11. Snapshots of a simulation with the Au(lll) surface Vads — 
10^s~^. Circles and diamonds denote fee and hep adsorption sites respec- 
tively. Filled and empty sites are represented by filled and empty figures 
respectively. The simulation time was 1.46 x 10~^ s. 



24 



This figure "Figurel.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/cond-mat/0106122vl 



This figure "Figure2.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/cond-mat/0106122vl 
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